RESONANT MESOMOLECULE FORMATION 
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Resonant formation of mcsomolecules in the mixture 
of hydrogen isotopes plays a key role in muon catalyzed 
fusion |l| . In the standard theory the rate of this pro- 
cess is calculated in the first order of perturbation the- 
ory, and perturbation potential V is calculated only in 
dipole approximation ||^, ^. The finite width of the 
resonance Q and the anharmonism of the potential |^ 
contribute essentially to the result. In the formation 
of the mesomolecule the most important perturbation 
is the shift of the deuteron ion in respect to the centre 
of mass [|[ : 



V{p,R = U{\p-P2R\)-U{p) 



(1) 



with potential U{p) for Z3X-molecules and mesomolec- 
ular complexes (MMCs). It determines the energy lev- 
els with different vibrational (i^) and rotational (K) 
quantum numbers and also calculates their wave func- 
tions. On the next stage the code calculates the ma- 
trix elements {vi, Ki, FlVli/f, Kf, S) where F is the 
spin of the original mesoatom and S is the spin of 
the mesomolecule. We need to know something like 
200 matrix elements and their dependence on the en- 
ergy. The code calculates also the position of all res- 
onances below and above the threshold. Next subrou- 
tine (SPEED) sums up all Breit-Wigner resonance con- 
tributions to the probability of mesomolecular complex 
formation A;^"*"^ and translates them from the centre- 



'\El,T)^ dE, 



Here U{p) is the ^Eg term of the DX mesomolecule of-mass frame to the lab. one 
{X = H,D,T); p is the distance between d^ and X+ 
ions, R is the distance between mesoatom and d^ in- 
side the mesomolecule, parameter /?2 (which depends on 
mass ratios) is /32=0.5-0.6. The potential U{p) is cal- 
culated in Born-Oppenheimer approximation (infinite 
mass of nuclei) and corrections to the finite nuclei mass 
are known to be small H- 



\E,)F{E,^El,T). 



Dipole interaction is only the first term of the multi- 
pole expansion for potential of Eq. (1). Meanwhile, the 
expansion parameter appears to be not a small num- 
ber. Indeed, this parameter is equal to the ratio of 
I32R to vibrational amplitude a = (Mfi)"^/^ which is 
/32i?(Mfi)i/2 ~ 0.2 - 0.35. For this reason the next 
(quadrupole) term of the multipole expansion changes 
formation rate significantly. 

In Ref. we developed a method which does not 
make use of the multipole expansion — the potential of 
Eq. (1) is taken into account exactly. We do not need 
the wave function of the mesomolecule, it is enough to 
know only its asymptotics. This is due to the fact that 
in the loosely bound mcsomolecules dyp (« = ,/=!; 
y = d,t) the mesoatom {yp)is is located in the average 
at large distances from the deuteron and hence one can 
neglect their interaction. Only exact binding energy 
I Ell I and constant Cdvfi which enters the asymptotics 
of the wave function are required for the calculation 
of formation rate. At the moment the binding energy 
I £11 1 is known with accuracy of the order of 1 meV. The 
most accurate values of the constant Cdy^i are obtained 
in Ref. pi. 



(2) 

Here the function F{Ec — > El,T), which takes 
into account the motion of gas molecules, leads to 
the Doppler broadening of resonances with the width 
Adq = {AE^^T/A)^^^ [|l|. In the spirit of the Bohr 
theory of compound nuclei, we consider the decay of 
mesomolecular complex as being independent of the 
channel by which it was formed. Therefore in order 
to obtain the rate of formation we have to multiply 
X^-^^{El,T) by the ratio Tr/Tt where Tt is the total 
and Tr is the partial reaction width (Tf for the fusion, 
^BD^ for the back decay, and so on). One has to 
average the values of A^ {E^ , T) over the spectrum of 
mesoatoms f{EL)- 

For the mesoatom spin F = we calculated the rates 
of resonant formation and fusion in reactions {tp.Y' +D2 
and {tpY' -\- DH at moderate gas density, when one can 
neglect rescattering of the MMC on the gas molecules 
[ p^ . In this case the decay of the MMC occurs from 
the state where it was formed. At small T the fusion 
rate is determined by underthreshold resonance with 
the energy Er — —U meV: 
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On the basis of the method developed in Q we elab- 
orated the computer code designed for the calculation 
of mesomolecule formation rates |l^. The subroutine 
MATEL of the code solves the Schrodinger equation 



that corresponds to the transition v = ^ 2, K = 
0^1. Here B = Noa{hc)^/me = 0.7085 • lO^^eV/c, 
f / = 0.65 meV and \{i\V\f)f = 0.60 • 10"* a.u. [|3|. 
According to Eq. (3), the rate A° does not depend 
on the temperature T. Account for other resonances 
leads to small variations of A° with temperature (see 
Fig. 1). This should be compared with the strong 
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Fi glire ll The rate of dtfi formation with subsequent fusion 
as a function of the temperature T at _F = exact (soUd 
curve) and in the approximation of zero width of resonances 
|l3| . Experimental point was measured at the density of equal 
to 1% of the density of liquid hydrogen. 



temperature dependence for the rate of dd^ formation 
(Fig. 2). Let us notice that, if one takes account of 
resonances above the threshold only (dashed line in 
Fig. 1), then the rate of formation restores its strong 
temperature dependence. Also let us point out that 
dipole approximation overestimates the matrix element 
of Eq. (1) more than 3 times 0, H). 

We plot on the Fig. 1 the single known experimental 
point. For the first time without any fitting parameter 
we reach the agreement with data on the level of 15% 
JTsf . The disagreement of previous models with data 
in the dtfi case was pointed out in a number of papers. 

The other limiting case when the collision width 
is much larger than all other widths takes place for 
the resonant formation of dd/x-mesomolecules in D2- 
gas [0. After formation of the mesomolecular com- 
plex in the state with rotational quantum number Kf 
it can be changed due to collisions to some new K'j:. 
It is enough time to reach the equilibrium Boltzman 
distribution in rotational quantum numbers K'j:. For 
this reason the back decay of the mesomolecular com- 
plex takes place independently of the value Kf with 
which it was formed. Rates of the ddfi formation with 
subsequent fusion for two values of the mesomolecule 
spin F = 1/2,3/2 are presented on Fig.2 ||l||. They 
are considerably less than for dtfi case. This fact 




Figure 2: Resonance rates of dd/i formation with subsequent 
fusion for mesoatoms {dfi)p with F = 1/2,3/2. Solid curves 
correspond to the theoretical calculations, dashed ones mark the 
error corridor ±4% for and ±3% for A^^^ . Experimental 
points were measured at the density equal to 5% of the density 
of liquid hydrogen. 



is a consequence of the sharp decrease of transition 

matrix elements with the energy of the bound state 

(e = -0.632 eV for dtfi and e = -1.966 eV for ddfj.). 

The rate A^^^(T) depends weakly on T owing to the 

strong Doppler broadening of the K = ^ 1 resonance 

with small positive resonance energy — 4.2 meV. 

1/2 

On the contrary, the rate (T) (which is mainly due 
to the far resonances with positive energy ) depends on 
the temperature rather strongly. 

At Fig. 2 we plot also recent experimental data of 
Ref. [1^ which were obtained for the non-equilibrium 
mixture of ortho- and para-deuterium in the ratio of 2 
to 1 (independent of T). In such a non-equilibrium gas 
due to the extra inelastic acceleration of mesoatoms 
on para-deuteuriums (which are excited by 7.4 meV), 
the mesoatom spectrum /{El) below « 70K should 
decrease with E^ considerably faster than the equilib- 
rium Maxwell spectrum fM{EL)- Meanwhile in our 
paper |l5j the fusion rate of Eq. (2) was averaged 
over fM{EL) spectrum. For this reason it makes sense 
to compare the results of Ref. with data only at 
temperatures T > 70K. The calculation (taking into 
account non-resonant formation of dd^ [[1^) based on 
purely theoretical values of all quantities, reproduces 
the data on the level of 10%. The accuracy of our cal- 
culations is limited mainly by insufficient accuracy of 
the estimation of fusion width F/ = 460(70) /isec~^ 
[ p^ . If one uses Fy as a single fitting parameter it is 
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Table 1: 





1 Cmatr \ 


P/, ^sec ^ 


lAE, 1993^^ |ig| 


0.62(13) 


337(40) 


lAE, 2001 |16| 


0.56(3) 


407(20) 


PNPI, 1998 [L5| 


1 


400(46) 


PNPL 2001 


1 


397(40) 



possible to achieve an excellent relative accuracy of 3% 
for A^/^(T) and 4% for \)''^{T). The best fit corre- 
sponds to P/ = 397(40) /isec~^. This value hes within 
the errors of theoretical estimate. 

Let us compare our results with calculations of Refs. 
[p^ , These calculations coincide with the data 

only if one introduces an additional fitting parameter 
\Cmatr\'^ = 0.56 which is used as factor to multiply all 
transition matrix elements squared. In other words, 
this approach overestimates by a factor of 1.8 (see 
Table 1). In our calculations this fitting parameter is 
absent. The discrepancy between the data and our 
calculations in the region of T < 7QK (up to 14% at 
T = 28.3 K) can be, in our opinion, attributed to the 
distortion of the spectrum /(E^) due to the mesoatom 
acceleration in non-equilibrium medium. This new ef- 
fect awaits for calculation. 

Thus the theory of the mesomolecule resonant forma- 
tion by Vesman mechanism developed in Refs. 

and improved significantly in Refs 0, |ll], |l^ de- 
scribes now the data with sufficient accuracy (see Figs. 
1,2). 
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